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Abstract. Ground-state behaviour of the frustrated quantum spin-i two- leg ladder 
with the Heisenberg intra-rung and Ising inter-rung interactions is examined in detail. 
The investigated model is transformed to the quantum Ising chain with composite 
spins in an effective transverse and longitudinal field by employing either the bond- 
state representation or the unitary transformation. It is shown that the ground 
state of the Heisenberg-Ising ladder can be descended from three exactly solvable 
models: the quantum Ising chain in a transverse field, the 'classical' Ising chain in 
a longitudinal field or the spin-chain model in a staggered longitudinal-transverse 
field. The last model serves in evidence of the staggered bond phase with alternating 
singlet and triplet bonds on the rungs of two-leg ladder, which appears at moderate 
values of the external magnetic field and consequently leads to a fractional plateau 
at a half of the saturation magnetization. The ground-state phase diagram totally 
consists of five ordered and one quantum paramagnetic phase, which are separated 
from each other either by the lines of discontinuous or continuous quantum phase 
transitions. The order parameters are exactly calculated for all five ordered phases 
and the quantum paramagnetic phase is characterized through different short-range 
spin-spin correlations. 
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1. Introduction 

Quantum spin models with competing interactions represent quite interesting and 
challenging topic for the current research [1, 2]. Such models show many unusual 
features in their ground-state properties and are very sensitive to the approximative 
schemes applied to them. An existence of exact solutions is therefore quite important, 
since they provide a rigorous information about the complex behaviour of frustrated 
models. However, exact results for such systems are still limited (see [3] and references 
cited therein). The simplest examples of rigorously solved quantum spin models are the 
models with the known dimerized [3-6] or even trimerized [7] ground states. Another 
example represents frustrated quantum spin antiferromagnets in high magnetic fields 
with the so-called localized-magnon eigenstates (see e.g. [8] for recent review). It should 
be nevertheless noticed that the aforementioned exact results are usually derived only 
under certain constraint laid on the interaction parameters and those models are not 
tractable quite generally. 

On the other hand, the exact solutions for some spin-^ quantum spin ladders 
with multispin interaction allows to find, besides the ground state, also thermodynamic 
properties quite rigorously [9,10]. Exactly solvable models, which admit the inclusion 
of frustration, can also be constructed from the Ising models where the decorations of 
quantum spins are included [11, 12]. The decoration-iteration procedure allows one to 
calculate exactly all thermodynamic properties of the decorated models when the exact 
solution for the corresponding Ising model is known. The distinctive feature of these 
models is that the quantum decorations are separated from each other, so that the 
Hamiltonian can be decomposed into the sum of commuting parts. The eigenstates of 
the total Hamiltonian are then simply factorized into a product of the eigenstates of its 
commuting parts. 

New findings of the exactly solvable frustrated quantum spin models with 
non-trivial (not simply factorizable) ground states are thus highly desirable. The 
present article deals with the frustrated quantum spin-| two- leg ladder with the 
Heisenberg intra-rung interaction and the Ising inter-rung interactions between nearest- 
neighbouring spins from either the same or different leg. Such a model can be 
regarded as an extension of the spin-^ Heisenberg-Ising bond alternating chain, which 
was proposed and rigorously solved by Lieb, Schultz and Mattis [13]. Alternatively, 
this model can also be viewed as the generalization of the exactly solved quantum 
compass ladder [14-16] when considering the Heisenberg rather than XZ-type intra-rung 
interaction and accounting for the additional frustrating Ising inter-rung interaction. 
Moreover, it is quite plausible to suspect that the exact results presented hereafter 
for the spin- 1 Heisenberg-Ising two-leg ladder may also bring insight into the relevant 
behaviour of the corresponding Heisenberg two-leg ladder, which represents a quite 
challenging and complex research problem in its own right [17,18]. The frustrated spin- 
| Heisenberg two-leg ladder have been extensively investigated by employing various 
independent numerical and analytical methods such as density-matrix renormalization 
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group [19-22], numerical diagonalization [22-26], series expansion [26], bosonization 
technique [27-29], strong- and weak-coupling analysis [30-34], valence-bond spin- wave 
theory [36], variational matrix-product approach [37,38], and bond mean-field theory 
[39-41]. Among the most interesting results obtained for this quantum spin chain, 
one could mention an existence of the columnar-dimer phase discussed in [32, 33] or a 
presence of the fractional plateau in the magnetization process examined in [22-25,34]. 

The theoretical investigation of two-leg ladder models is motivated not only from 
the academic point of view, but also from the experimental viewpoint, because the 
two-leg ladder motif captures the magnetic structure of a certain class of real quasi- 
one-dimensional magnetic materials. The most widespread family of two-leg ladder 
compounds form cuprates, in which one finds both experimental representatives with 
the dominating intra-rung interaction like SrCu203 [42,43], Cu2(C5Hi 2 N 2 )2Cl4 [44-47], 
(C5Hi 2 N) 2 CuBr4 [48], (5IAP) 2 CuBr 4 [49] as well as, the magnetic compounds with the 
dominating intra-leg interaction such as KCuCl 3 [50,51], TlCuCl 3 [52-54], NH 4 CuCl 3 , 
KCuBr 3 . Understanding the low-temperature magnetism of two- leg ladder models also 
turns out to be crucial for an explanation of the mechanism, which is responsible for the 
high-temperature superconductivity of cuprates [55]. Besides the cuprates, another 
experimental representatives of the two-leg ladder compounds represent vanadates 
(VO) 2 P 2 7 [56], CaV 2 5 and MgV 2 5 [57], as well as, the polyradical BIP-BNO [58,59]. 

The outline of the paper is as follows. In section 2 the model is defined and the 
pseudospin representation is considered. The ground-state properties of the model with 
and without external field are studied in section 3. The most important findings are 
summarized in section 4. 

2. Heisenberg-Ising two-leg ladder 

Let us define the spin-| Heisenberg-Ising ladder through the following Hamiltonian: 



N 



H 



y] \Jl( S l,i ■ S 2,i)a + <MSl,tSl,i+l + S 2,i S 2,i+l) 



8=1 



+ J-i{ S 2,i S l,i+l + S l,i S 2,i+l) ~ H S l,i + S 2,i) 



(1) 




2,N 



l.N 



Figure 1. Quantum spin-i Heisenberg-Ising two-leg ladder. Thick (thin) lines denote 
the Heisenberg (Ising) bonds. 
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|4_> = -^(ItMtw) " |<k*ki», lft,+> = -sCltMtw) + Ui,4a,i». (2) 



where (s^j -s 2i j)a = s i,j s 2,i + s i i s 2 i + AsfjS 2i , s h are three spatial projections of spin-| 
operator, the first index denotes the number of leg, the second enumerates the site, J\ 
is the XXZ Heisenberg intra-rung interaction between nearest-neighbour spins from 
the same rung, J 2 is the Ising intra-leg interaction between nearest-neighbour spins 
from the same leg, J3 is the crossing (diagonal) Ising inter-rung interaction between 
next-nearest-neighbour spins from different rungs, h is the external magnetic field. We 
also imply the periodic boundary conditions s 1A r +1 = s 1;1 and s 2)A r+i = s 2 ,i along legs. 
The coupling constants J 2 and J3 can be interchanged by renumbering of the sites, as 
well as their signs can be simultaneously reverted by spin rotations. Therefore, the 
Hamiltonians H[J 2 , J3), H{J 3 , J 2 ) and H(— J 2 , — J3) in zero field have equal eigenvalues 
and the corresponding models are thermodynamically equivalent. 

It can be checked that ^-projection of total spin on a rung Sf = sf j + s 2 i 
commutes with the total Hamiltonian [Sf, H] = and hence, it represents a conserved 
quantity. For further convenience, it is therefore advisable to take advantage of the 
bond representation, which has been originally suggested by Lieb, Schultz and Mattis 
for the Heisenberg-Ising chain [13]. Let us introduce the bond-state basis consisting of 
four state vectors: 

Ko> = -^={\U,ih,i) ~ ITmW), \4>{,o) = -^(| Wt2,i) + lti,iW)> 

1 1 

-^(Iti/la,*) - Ui,4 2 ,i», = ^ 

These states are the eigenstates of the Heisenberg coupling between two spins located 
on the ith rung, i.e. |0o,o) ^ s the singlet-bond state, \(j>\ ), \(j>i ±) are triplet states. 
Following [13] two subspaces can be singled out in the bond space: i) \<f> l 0Q ), \4>\o) and 
ii) + ), _). The first and second subspaces corresponds to {Sf) 2 = and {Sf ) 2 = 1 
respectively, and the Hamiltonian can be diagonalized separately in each subspace. We 
can introduce the index of subspace at i-th site: n, = 0(1) if the given state is in 
the subspace with {Sf) 2 = {{Sf) 2 = 1) spanned by states |0o,o)> I0i,o) G0i,-)> I0i,+))- 
According to [13], let us also call the bond states as purity (impurity) states if = 0(1). 
If one makes pseudospin notations for the states |0oo) = H)o> I0io) = lt)o> the action 
of the spin- 1 operators in the subspace = can be expressed in terms of new raising 
and lowering operators: 

1 1 

4,i = -^(af + at) = Si, 4,i = + Oi) = sf, 

, 2 + A _ A 
(si,i • s 2ji ) A = aja t — = s { - — . (3) 

The operators af, aj satisfy the Pauli algebra {{ai,af} = 0, {ai,ai} = {af,af} = 0, 
[cij, aj] = [a^ at] = [af , at] = for i / j), and half of their sum can be identified as a 
new pseudospin operator sf. 

Analogously, one can consider the pseudospin representation of the subspace = 1 
(|0i _} = |4) 1, 101+) = |t)i) and find the action of spin operators in it as follows: 

4,i = + <h) = sf, s z %i = ^(a+ + af) = sf, (s M ■ s 2ji ) A = j- (4) 
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Combining (3), (4) we find the general expressions for the pseudospin representation of 
these operators, which are valid in both subspaces: 

4,i = (2ni-l)sf, s z %t = sf, (si, i -s 2)i )A = (l-n i )a? + j(2ni-l). (5) 
The effective Hamiltonian can be rewritten in terms of new operators as follows: 



H 



£ { Ji(l - ni)5f - 2hn i r i + - J ^(2n i - 1) 

i=l 



+ [4J 2 n^ +1 - 2(J 2 - J 3 ){m + n m - 1)] S?S? +1 }. (6) 

It is noteworthy that equation (5) can be considered as some nonlinear spin 
transformation from sf i} s%i to new operators s'"^ s'^, where e.g. s'^ = (rij — ~), 
s?2i — sf. The transformation can be generated by the following unitary operator: 

N 

U = Y[exp + ex P exp ( _i I S 2,i) ex P ( i7r 4<) " ( 7 ) 

i=l 

The crucial point is the second factor, which introduces the nonlinearity. Other terms 
are the rotation in spin space, they are used to adjust the transformed operators to the 
form of (5). Finally, the spin operators are transformed as follows: 

Ust,U + = aft, UslU + = Usi t U + = 2^, 

E^IT* = 2^^, = - 2s iA,i, UslJJ + = si,. (8) 

The transformed Hamiltonian has the form of the quantum Ising chain with composite 
spins in an effective longitudinal and transverse magnetic field: 

N J J A 

UHU+ = { -j (1 - 2*y 4, + ~~2~ s i,i + [M^iA+i + 1) 
1=1 

+ 2J 8 (a? fi+1 + st,)] 4,4,.: - Ml + ( 9 ) 

The straightforward correspondence sf j = n.j — | , s 2 j = sf leads to the equivalence 
between the Hamiltonians (6) and (9). It should be noted that such kind of a 
representation remains valid also in case of the asymmetric ladder having both diagonal 
(crossing) Ising interactions different from each other. 



3. Ground state of the Heisenberg-Ising two-leg ladder 

The Hamiltonian (6) can also be rewritten in the following more symmetric form: 

N j ^ 

H= £ [ Jx(l - m)st - 2hn l ~s x i + ^(2n, - 1) 
i=i 

+ [2(J 2 + J 3 )n ini+1 + 2(J 2 - J 3 )(l - m)(l - n i+l )\ (10) 

This form serves in evidence that the effective model splits at zth site into two 
independent chains provided that two neighbouring bonds are being in different 
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subspaces, i.e. 7^ n i+ i. The uniform Hamiltonians, when all = or all = 1, 
read: 

N 



£ [ Jl <5 " t) + 2( J2 " J3 ^ S ~M ' 

1 = [ 2 (^ + - 2/w? + ^] . (12) 



JV 
i=l 



If all bonds are in purity states (rij = 0), one obtains the effective Hamiltonian of 
the Ising chain in the transverse field that is exactly solvable within Jordan- Wigner 
fermionization [60,61]. If all bonds are in impurity states (n, = 1), one comes to the 
Ising chain in the longitudinal field that is solvable by the transfer-matrix method (see 
e.g. [62]). Accordingly, the ground-state energy of the model in = subspace is given 
by [61]: 

eS Um Sp. = J^i J -- J \ (v rf) - M, (13) 

TV — s-oo iV 7T 4 

where 7 = ^ j^~j 3 and E(ac) = / d9^/l — k 2 sin 2 9 is the complete elliptic integral of 
the second kind. 

The ground-state energy of the model in rij = 1 subspace is given by the ground- 
state energy of the effective Ising chain in the longitudinal field: 

6 °-nZo N ~ \ - *±*, if \h\ < (J 2 + J S ). U4j 

The upper (lower) case corresponds to the ferromagnetically (antiferromagnetically) 
ordered state being the ground state at strong (weak) enough magnetic fields. 

3.1. Ground-state phase diagram in a zero field 

If the external field vanishes (h = 0), it is sufficient to show that the inequality 
Eq(Ni) + Eq(N 2 ) > Eq(Ni + N 2 ) holds for any finite transverse Ising chain with free 
ends in order to prove that the ground state always corresponds to a uniform bond 
configuration. Here E®(N) denotes the ground state energy of the Hamiltonian: 

N-l N . 

H°(N) = 2( J 2 - J 3 ) s >7 + i + { S i - t) • ( 15 ) 



i=l i=l 

Indeed, two independent chains of size Ni and N 2 can be represented by the following 
Hamiltonian: 

iVi-l JVi+JVa-l N!+N 2 



A 

7 



H°{N U N 2 ) = 2( J 2 - J 3 ) £ SfS^ + 2( J 2 - J S ) £ S?S^ + ( s 1 

i=l i=JVi+l i=l 

= + iV 2 ) - 2(J 2 - J 3 )~s x Nl ~s x Nl+1 . (16) 

If Eq(N) and |^) are the lowest eigenvalue and eigenstate of H°(N), then E®(Nx, N 2 ) = 
Ef&Nj + E${N 2 ) and \i(>o uNa ) = itf 1 ) W a ) are the lowest eigenvalue and eigenstate of 
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H ( N\ , N 2 ) . Now, it is straightforward to show that 
£ °(Aq) + E° (N 2 ) = (^IH^N, + N 2 )\^) 

- 2(J 2 - J 3 )(^ 1 |^ 1 K 1 )« 2 |^ 1+1 K 2 ) > E^N, + N 2 ). (17) 
Here, we have used that (V'o' 1 \s% 1 li^o 1 ) = for any finite chain [13,61], and that the 
lowest mean value of the operator H°(Ni + N 2 ) is achieved in its ground state. 

It is easy to show by straightforward calculation that the same property is valid for 
the Ising chain with free ends in zero longitudinal field, i.e. 

J5£(JVi + N 2 ) = E^N,) + E*(N 2 ) + l A±^A < E ^ Nl) + jgg^). (18) 

Now, let us prove that the ground state of the whole model may correspond 
only to one of uniform bond configurations. It can be readily understood from (10) 
that the effective Hamiltonian does not contain an interaction between spins from 
two neighbouring bonds if they are in different subspaces [63]. It means that the 
effective model splits into two independent parts at each boundary ('domain wall') 
between the purity and impurity states. Thus, the Heisenberg-Ising ladder for any 
given configuration of bonds can be considered as a set of independent chains of two 
kinds and of different sizes. Then, the ground-state energy of any randomly chosen bond 
configuration will be as follows: 

E = £°(jV!) + E° (N 2 ) + ... + El(Mi) + El(M 2 ) + ... 

> Eg(JV! + N 2 + ...) + E^{M 1 + M 2 + . . .) (19) 

It is quite evident from equations (17) and (18) that the one uniform configuration, 
which corresponds to the state with lower energy than the other one, must be according 
to inequality (19) the lowest-energy state (i.e. ground state). Therefore, the model 
may show in the ground state the first-order quantum phase transition when the lowest 
eigenenergies in both subspaces becomes equal e\ = e$. Besides, there also may appear 
the more striking second-order (continuous) quantum phase transition at | J\\ = \ J 2 — J$\ 
in the ground state, which is inherent to the pure quantum spin chain [61]. From this 
perspective, the Heisenberg-Ising ladder can show a variety of quantum phase transitions 
in its phase diagram. 

The ground-state phase diagram of the spin-| Heisenberg-Ising ladder in a zero 
magnetic field is shown in figure 2. It reflects the symmetry of the model, i.e. it is 
invariant under the exchange of the Ising interactions J 2 and J3, as well as, under the 
simultaneous change of the signs of J 2 and J 3 . Altogether, five different ground states 
can be recognised in figure 2: 

(i) Quantum paramagnetic (QPM) state for e° < ej, and \J\\ > \ J 2 — J 3 |: the 
equivalent transverse Ising chain (11) is in the gapped disordered state with no 
spontaneous magnetization (sf) = and non-zero magnetization (sf) =/ induced 
by the effective transverse field. For the initial Heisenberg-Ising ladder it means that 
the rung singlet dimers are dominating on Heisenberg bonds in the ground state. 
Although each bond may possess the antiferromagnetic order, the interaction along 
legs demolishes it and leads to the disordered quantum paramagnetic state. 




Figure 2. Ground-state phase diagrams for the spin-^ Hciscnbcrg-Ising two- leg ladder 
in an absence of the external magnetic field (h = 0). Broken (solid) lines denote the 
lines of first-order (second-order) quantum phase transitions, dotted lines denote the 
rung singlct-dimcr state. 

(a) Ji = 1, A = O.O(red), l.O(blue), 1.5(green) (from thick to thin lines); 

(b) J% = 1, A = 1. Thick solid and broken (blue) lines are the ground-state boundaries 
of the Heisenberg-Ising ladder, the thin line with dots shows the ground-state boundary 
of the Heisenbcrg ladder [26] and straight broken lines the ground-state boundaries of 
the Ising ladder. 



(ii) Stripe Leg (SL) state for e[] < ej, and | J\\ < J3 — J2: the equivalent transverse 
Ising chain exhibits the spontaneous ferromagnetic ordering with (sf) 7^ 0. Due to 
relations (3) and (4), one obtains for the Heisenberg-Ising ladder (sf j) = (sf = 
~( s 2i) = ~( s 2i+i) 7^ 0- This result is taken to mean that the Heisenberg-Ising 
ladder shows a ferromagnetic order along legs and antiferromagnetic order along 
rungs, i.e. the magnetizations of chains are opposite. The staggered magnetization 
as the relevant order parameter in this phase is non-zero and it exhibits evident 
quantum reduction of the magnetization given by: m z SL = ^ ^2iLii( s i i) ~ ( s 2 i)) = 
|[1 — J\j(Ji — Jz) 2 ^- Here we used the result for the spontaneous magnetization 
of transverse Ising chain [61]. 

(iii) Neel state for e[j < ej, and |Ji| < J2 — J3'- the effective transverse Ising chain 
shows the spontaneous antiferromagnetic order with (sf) = (—l) t m x ^ 0. For the 
Heisenberg-Ising ladder one consequently obtains (sf j) = — (sf = — (sfi) = 
(s2i+i) 7^ 0- Hence, it follows that the nearest-neighbour spins both along legs as 
well as rungs exhibit predominantly antiferromagnetic ordering. The dependence 
of staggered magnetization as the relevant order parameter is quite analogous to 
the previous case m z AF = ± E J =i(- 1 ) i (( s ti) ~ ( s 2,i)) = It 1 ~ J i/( J 2 ~ Js) 2 }*- 

(iv) Stripe Rung (SR) state for < e[j, J2 > 0, J3 > 0: the model shows classical 
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ordering in this phase with the antiferromagnetically ordered nearest-neighbour 
spins along legs and the ferromagnetically ordered nearest-neighbour spins along 
rungs. 

(v) Ferromagnetic (FM) state for < e{], J 2 < 0, J3 < 0: the ground state 
corresponds to the ideal fully polarized ferromagnetic spin state. 

The results displayed in figure 2(a) demonstrate that the Heisenberg-Ising ladder 
is in the disordered QPM phase whenever a relative strength of both Ising interactions 
J 2 and J3 is sufficiently small compared to the Heisenberg intra-rung interaction J\. It 
is noteworthy, moreover, that the QPM phase reduces to a set of fully non-correlated 
singlet dimers placed on all rungs (the so-called rung singlet-dimer state) along the 
special line J 2 = J3 up to < ej, which is depicted in figure 2 by dotted lines. Under 
this special condition, the intra-rung spin-spin correlation represents the only non-zero 
pair correlation function and all the other short-ranged spin-spin correlations vanish 
and/or change their sign across the special line J2 = J3. It should be remarked that the 
completely identical ground state can also be found in the symmetric Heisenberg two-leg 
ladder with J 2 = J3 (see figure 2(b)). To compare with, the rung singlet-dimer state 
is being the exact ground state of the symmetric Heisenberg-Ising ladder with A = 1 
for J 1 > J 2 , while the symmetric Heisenberg ladder displays this simple factorizable 
ground state for J\ > 1.4015 J2 [35,36] (this horizontal line is for clarity not shown in 
figure 2(b) as it exactly coincides with the ground-state boundary of the Heisenberg- 
Ising ladder extended over larger parameter space). It should be stressed, however, 
that the short-range spin-spin correlations become non-zero in QPM whenever J 2 7^ J3 
even if this phase still preserves its disordered nature with the prevailing character of the 
rung singlet-dimer state. To support this statement, the zero-temperature dependencies 
of the order parameters and the nearest-neighbour spin-spin correlation along legs are 
plotted in figure 3. The relevant order parameters evidently disappear in QPM, whereas 
the nearest-neighbour correlation function changes its sign when passing through the 
special line J 2 = J3 of the rung singlet-dimer state. 

Note furthermore that the Heisenberg-Ising ladder undergoes the second-order 
quantum phase transition from the disordered QPM phase to the spontaneously long- 
range ordered Neel or SL phase, which predominantly appear in the parameter region 
where one from both Ising couplings J 2 and J3 is being antiferromagnetic and the other 
one is ferromagnetic. The significant quantum reduction of staggered order parameters 
implies an obvious quantum character of both Neel as well as SL phases in which the 
antiferromagnetic correlation on rungs still dominates (see figure 3). Thus, the main 
difference between both the quantum long-range ordered phases emerges in the character 
nearest-neighbour correlation along legs, which is ferromagnetic in the SL phase but 
antiferromagnetic in the Neel phase. The continuous vanishing of the order parameters 
depicted in figure 3(a) provides a direct evidence of the second-order quantum phase 
transition between the disordered QPM phase and the ordered SL (or Neel) phase, 
which is also accompanied with the weak-singular behaviour of the nearest-neighbour 
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Figure 3. Zero-temperature variations of: (a) the order parameters; (b) the nearest- 
neighbour spin-spin correlation along legs; as a function of the Ising intra-leg interaction 
Ja for Ji = 1, A = 1 and three different values of the crossing Ising interaction 
J 3 = 0,0.6,1-2. 



correlation function visualized in figure 3(b) by black dots. Here it should be noted 
that the zz correlation function of the Heisenberg-Ising ladder can be easily derived 
from the result of the xx correlation function of the transverse Ising chain calculated 
in [61]. Last but not least, the strong enough Ising interactions J 2 and J 3 may break 
the antiferromagnetic correlation along rungs and lead to a presence of the fully ordered 
ferromagnetic rung states FM or SR depending on whether both Ising interactions are 
ferromagnetic or antiferromagnetic, respectively. It should be noticed that this change 
is accompanied with a discontinuous (first-order) quantum phase transition on behalf 
of a classical character of both FM and SR phases, which can be also clearly seen in 
figure 3 from an abrupt change of the order parameters as well as the nearest-neighbour 
correlation function. 

Finally, let us also provide a more detailed comparison between the ground- 
state phase diagrams of the Heisenberg, Heisenberg-Ising and Ising two-leg ladders 
all depicted in figure 2(b). One can notice that the displayed ground-state phase 
diagrams of the Heisenberg and Heisenberg-Ising two-leg ladders have several similar 
features. The phase boundary between the classically ordered SR phase and three 
quantum phases (QPM, SL, Neel) of the Heisenberg-Ising ladder follows quite closely 
the first-order phase boundary between the Haldane-type phase and the dimerized 
phase of the pure Heisenberg ladder obtained using the series expansion and exact 
diagonalization [26]. Of course, the most fundamental difference lies in the character 
of SR and Haldane phases, because the former phase exhibits a classical long-range 
order contrary to a more spectacular topological order of the pure quantum Haldane- 
type phase with a non-zero string order parameter [31] even though the ferromagnetic 
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intra-rung correlation is common for both phases. The difference between the ground 
states of the Heisenberg-Ising and pure Heisenberg two-leg ladder becomes much less 
pronounced in the parameter region with the strong intra-rung interaction Ji, which 
evokes the quantum phases in the ground state of both these models. In the case of a 
sufficiently strong frustration J 2 ~ J3, the ground state of the Heisenberg-Ising ladder 
forms the disordered QPM phase, whereas the antiferromagnetic or ferromagnetic (Neel 
or SL) long-range order emerges along the legs if the diagonal coupling J 3 is much weaker 
or stronger than the intra-leg interaction J 2 . Note that the quantum Neel and SL phases 
have several common features (e.g. predominant antiferromagnetic correlations along 
rungs) with the disordered QPM phase from which they evolve when crossing continuous 
quantum phase transitions given by the set of equations: J 3 = J 2 _ J\ (for J 3 < J 2 ) 
and J3 = J 2 + Ji (for J3 > J 2 ). It cannot be definitely ruled out whether or not these 
quantum ordered phases may become the ground state of the pure Heisenberg ladder, 
because they are also predicted by the bond-mean-field approximation [39-41] but have 
not been found by the most of the numerical methods. Further numerical investigations 
of the Heisenberg ladder are therefore needed to clarify this unresolved issue. 



3.2. Ground-state phase diagram in a non-zero field 

When the external field h is switched on, the inequality for the ground states of the 
Ising chain in the longitudinal field (18) is generally valid only for J\ + J 2 < or 
J\ + J2 > 2\h\. Hence, it is necessary to modify the procedure of finding the ground 
states inside the region where the relation (18) is broken. However, one may use with 
a success the method suggested by Shastry and Sutherland [25,64] in order to find the 
ground states inside this parameter region. Let us represent our effective Hamiltonian 
(10) in the form: 

N 

h = y~]Hi )i+ i, 

i=l 

H i>i+1 = - { -Ml - niYst - 2hntsf + \{2n, - 1)} 



+ [2(J 2 + J 3 )mn i+1 + 2(J 2 - J 3 )(l - m){l - n i+x )\ S?S? +1 . (20) 

Then, one can employ the variational principle implying that E > X^i-^o(M + 1)j 
where E (l, I + 1) is the lowest eigenenergy of H^i +1 . 

Looking for the eigenenergies of i^z+i it is enough to find the lowest eigenstate of 
each bond configuration: 

j< > (4 -n - _1. / 72 _l ( Z _ JA2 _ JiA- 

Hi = n i+1 = 1, 



Hi = n i+1 = 0, Eq ,0 (i,i + 1) = -|a/Ji + (J 2 - J3) 




FV(ii + H-l ^r + ^-H if|/i|>(J 2 + J 3 )(FM), 

MMHiJ-i iA,^ if|/i|<(J2 + Js)(AF). 
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• m = o, m+i = 1{th = 1, m+i = 0), + l) = E 1 '°(i,i + 1) = -f - M 

The phase corresponding to the alternating bond configuration = 0, = 1 will 
be hereafter referred to as the staggered bond (SB) phase. It should be mentioned 
that there does not exist in the SB phase any correlations between spins from different 
rungs and the overall energy comes from the intra-rung spin-spin interactions and the 
Zeeman's energy of the fully polarized rungs [63]. It can be easily seen from a comparison 
of E^ii, 2 + 1) and E ' (i, i + 1) that the eigenenergy of the SB phase Eq' (i, i + 1) has 
always lower energy than the lowest eigenenergy of the fully polarized state E ' (i, i + 
inside the stripe: 

j, + J ,-!^< w <j, + j, + W + *L. (21) 

Note that the ferromagnetic ordering is preferred for external fields above this stripe 
\h\ > J 2 + ^3 + Jl ( 1 2 +A - ) , while the antiferromagnetic ordering becomes the lowest-energy 

7 3 - -^T 



state below this stripe \h\ < J 2 + 



If one compares the respective eigenenergies of the staggered bond phase Eq' 1 ^, 2+1) 
and the uniform purity phase E '°(i,i + 1) (n, = n i+1 = 0), one gets another condition 
implying that Eq' (i, i + becomes lower than EQ ,0 (i, i + 1) only if: 



\h\ > y/j! + ( J 2 - J 3 ) 2 - Jl(1 2 A) . (22) 

This means that the SB phase with the overall energy Eq' 1 = —N(Ji + 2|/i|)/4 becomes 
the ground state inside the region confined by conditions (21) and (22). The lowest 
field, which makes the SB phase favourable, can be also found as: 

Similarly, one may also find the condition under which two eigenenergies 
corresponding to the uniform impurity configuration E ' (rii = rij + i = 1) become lower 
than the respective eigenenergy of the uniform purity configuration E ' (rii — n i+i — 0). 
The ferromagnetic state of the uniform impurity configuration becomes lower if the 
external field exceeds the boundary value: 



i i - 2 



J 2 + J 3 + J 1 A + ,/J 1 2 + (j 2 _ j 3 ; 



(24) 



whereas the condition for the antiferromagnetic state of the uniform impurity 
configuration is independent of the external field: 



JxA + y/J* + (J 2 - J 3 ) 2 < (J 2 + J 3 ). (25) 

It is worthy of notice that it is not possible to find the ground state outside the 
boundaries (22), (24), (25) using the variational principle. However, it is shown in the 
appendix that the bond configuration, which corresponds to the ground state, cannot 
exceed period two. Therefore, one has to search for the ground state just among the 
states that correspond to the following bond configurations: all rii = 0; all rii — 1; 
ri2i-\ = 0, n 2 i = 1 [ri2i-i = 1, ri2i = 0). In this respect, two phases are possible inside 
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Figure 4. Ground-state phase diagrams for the Heisenberg-Ising two-leg ladder in 
the J2 — J3 plane for Ji = 1, A = 1 and two different values of the external field: (a) 
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= 1.0; (b) h = 1.1. 
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Figure 5. Ground-state phase diagrams for the Heisenberg-Ising two-leg ladder in the 
J2 — h plane for Ji = 1, A = 1 and two different values of the crossing interaction: (a) 
J 3 = 0; (b) J 3 = 1. 



the stripe given by (21). The ground state of the Heisenberg-Ising ladder forms the SB 
phase if e[j' < e[j: 

> 2(J 1 + IW3l) E(yr3?) _ Ji(l-A)_ ( 

7T 2 

Several ground-state phase diagrams are plotted in figures 4-6 for a non-zero 
magnetic field. The most interesting feature stemming from these phase diagrams is 
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that the SB phase may become the ground state for the magnetic field h > Jl ( 1 + A ) ; 
which is sufficiently strong to break the rung singlet-dimer state. It can be observed 
from figure 4 that the SB phase indeed evolves along the line of the rung singlet-dimer 
state and replaces the QPM phase in the ground-state phase diagram. The external 
magnetic field may thus cause an appearance of another peculiar quantum SB phase 
with the translationally broken symmetry, i.e. the alternating singlet and fully polarized 
triplet bonds on the rungs of the two-leg ladder. Hence, it follows that the SB phase 
emerges at moderate values of the external magnetic field and it consequently leads 
to a presence of the intermediate magnetization plateau at a half of the saturation 
magnetization. It is quite apparent from figures 5(a), 6(a) that the Heisenberg-Ising 
ladder without the frustrating diagonal Ising interaction J 3 = exhibits this striking 
magnetization plateau just for the particular case of the antiferromagnetic Ising intra-leg 
interaction J2 > 0. On the other hand, the fractional magnetization plateau inherent to 
a presence of the SB phase is substantially stabilized by the spin frustration triggered by 
the non-zero diagonal Ising interaction J3 7^ and hence, the plateau region may even 
extend over a relatively narrow region of the ferromagnetic Ising intra-leg interaction 
J 2 < as well (see figure 5(b)). 

It should be noted that the same mechanism for an appearance of the magnetization 
plateau has also been predicted for the spin-| Heisenberg two-leg ladder by making use 
of exact numerical diagonalization and DMRG methods [22-25,34]. Let us therefore 
conclude our study by comparing the respective ground-state phase diagrams of the 
Heisenberg-Ising and Heisenberg two-leg ladders in a presence of the external magnetic 
field displayed in figure 6. According to this plot, both models give essentially the 
same magnetization process in the strong-coupling limit of the Heisenberg intra-rung 
interaction J\ <; 1.5, where two subsequent field-induced transitions in the following 
order QPM-SB-FM can be observed and consequently, the magnetization exhibits two 
abrupt jumps at the relevant transition fields. On the other hand, the fundamental 
differences can be detected in the relevant magnetization process of the Heisenberg- 
Ising and Heisenberg ladders in the weak-coupling limit of the intra-rung interaction 
J 1. At low fields, the quantum Haldane-like phase constitutes the ground state of the 
pure Heisenberg ladder in contrast to the classical SR phase, which is being the low-field 
ground state of the Heisenberg-Ising ladder. In addition, the Heisenberg-Ising ladder still 
exhibits a magnetization plateau corresponding to a gapped SB phase at intermediate 
values of the magnetic field in this parameter space, while a continuous increase of the 
magnetization can be observed in the pure Heisenberg ladder due to a presence of the 
gapless Luttinger-liquid phase at moderate fields. Finally, it is also noteworthy that the 
saturation fields towards FM phase are also quite different for the Heisenberg-Ising and 
Heisenberg ladders in the weak-coupling limit of the Heisenberg intra-rung interaction. 




Figure 6. Ground-state phase diagram for the Heisenberg-Ising two- leg ladder in the 
J\ — h plane for J2 = 1, A = 1 and two different values of the crossing interaction: (a) 
J3 = 0; (b) J 3 = 1. For a comparison, figure 6(b) depicts by red lines the ground-state 
phase diagram of the pure Heisenberg two-leg ladder adapted from [22]. 



4. Conclusions 

The frustrated Heisenberg-Ising two-leg ladder was considered within the rigorous 
approach using that ^-projection of total spin on a rung is a conserved quantity. 
By means of the pseudospin representation of bond states, we have proved the exact 
mapping correspondence between the investigated model and some generalized spin-^ 
quantum Ising chain with the composite spins in an effective longitudinal and transverse 
field. We have also found the unitary transformation which reproduces the analogous 
mapping between the models. While the bond-state representation has more transparent 
physical interpretation, the unitary transformation gives the complete relation between 
spin operators of both models, and it might be useful for searching quantum spin ladders 
which admit exact solutions. 

It has been shown that the ground state of the model under investigation must 
correspond to a regular bond configuration not exceeding the period two. Hence, the 
true ground state will be the lowest eigenstate of either the transverse Ising chain, 
the Ising chain in the longitudinal field, or the non-interacting spin-chain model in the 
staggered longitudinal-transverse field. 

The most interesting results to emerge from the present study are closely related 
to an extraordinary diversity of the constructed ground-state phase diagrams and the 
quantum phase transitions between different ground-state phases. In an absence of the 
external magnetic field, the ground-state phase diagram constitute four different ordered 
phases and one quantum paramagnetic phase. The disordered phase was characterized 
through short-range spin correlations, which indicate a dominating character of the 
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rung singlet- dimer state in this phase. On the other hand, the order parameters have 
been exactly calculated for all the four ordered phases, two of them having classical 
character and two purely quantum character as evidenced by the quantum reduction of 
the staggered magnetization in the latter two phases. Last but not least, it has been 
demonstrated that the external magnetic field of a moderate strength may cause an 
appearance of the peculiar SB phase with alternating character of singlet and triplet 
bonds on the rungs of two-leg ladder. This latter finding is consistent with a presence 
of the fractional magnetization plateau at a half of the saturation magnetization in the 
relevant magnetization process. 
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Appendix A. Ground-state bond configuration in magnetic field 

In general each state of the model can be represented as an array of alternating purity 
and impurity non-interacting clusters of different length. One can formally write the 
lowest energy of some configuration as 

E = E (N U Mi) + E (N 2 , M 2 ) + . . . + E (N L , M L ), (A.l) 

where E (Ni, M { ) = E°(Ni) + Eq (Mi) is the lowest energy of the complex of two 
independent spin chains which correspond to the purity and impurity bond states, 
N 1 + M 1 +N 2 + M 2 + . ..N L + M L = N. If e Q (N u Mi) = E (Ni, M i )/(N i +M i ) corresponds 
to the lowest energy per spin among other clusters, it is evident that the ground state 
configuration is alternating purity and impurity clusters of length iVj and Mj. 

Let us consider at first the case when h is restricted by condition (21). If the 
configuration contains the cluster of more than 2 impurity bonds, its energy can be 
lowered by adding pure bond in-between. If condition (21) is valid, such a configuration 
can achieve a lower energy. Using this procedure we can reduce the number of sites 
in impurity clusters to 2, i.e. only E (N, 1) or E (N,2) can correspond to the ground 
state. 

Define the energy per spin for each configuration as: 

^D ^^f 1 , (A.2) 

eo(N , 2)= mq±m. (A .3) 

If e (l, 1) < e (2, 1) and so on (e (iV, 1) < e (iV + 1, 1)) the ground state corresponds 
to the staggered bond configuration. Suppose that for some N 

Ae (iV-l)=e (iV,l)-e (iV-l,l) 
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Figure Al. A£g(2V) = E$(N + 1) - £§(iV) as a function of system size N for the 
transverse Ising chain with free ends (15): J\ = 1, A = 0, | J2 — J3I = 0, 0.2, 0.4, . . . , 4 
from top to bottom. 



N(N 



T (NE° (N) - (N + 1)E° (N - 1) - ^(1)) < 0. (A.4) 



Consequently, it is clear that 

Aeo(iV) = (Ar + 1) 1 (Ar + 2) ((iV + 1)E° Q (N + 1)-(N + 2)E° Q (N) - J5*(l)) < 0, (A.5) 

if Eq(N + 1) — 2Eq(N) + Eq(N — 1) < 0, i.e. the ground state energy of finite transverse 
Ising chain is a convex function of N. The last condition is ensured by non-increasing 
AEq(N) = E$(N+1)-E$(N) (AEq(N) > AE$(N+1)) that can be shown by numerical 
calculations for finite chains (see figure Al) using the numerical approach [65]. For 
J 2 = J3 we get the free spin Hamiltonian and E®(N) depends linearly on N. Conditions 
(A. 4), (A. 5) together mean that eo(N, 1) may have only one extremum and it is the 
maximum. Therefore, eo(A^, 1) may take the minimal value only in two limiting cases 
N = 1, N — > 00. The same result can be analogously obtained for eo(A r , 2). Now we 
can compare the energies of (N, 1) and (N,2) configurations. For \h\ > Jl \ j2 

e (l, 1) - e (l, 2) = i (- Jl(1 2 +A) - ( J 2 + J 3 ) + \h\\ < 0, (A.6) 

due to the upper boundary of condition (21). To summarize, the ground state 
configuration in the region confined by (21) can correspond to either staggered bond or 
uniform purity configuration. 

Let us consider the ground state for \h\ > J2 + J3 + . Similarly to the previous 

case we can prove that e (A^, M) is a convex function of and M if the ground state 
energies of the uniform chains have the following properties: 

E°(N + 1) - 2E° (N) + E°(N - 1) < 0, 

Eq(N + 1) - 2EKN) + E^N - 1) < 0. (A.7) 
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One can find that E%(N+1)-E%(N) = ^f^, and E^(N+1)-2E^(N)+E^(N-1) = 0. 
That is why the minimal value of the ground state can be achieved in one of three cases: 
staggered bond, purity and impurity configuration. The staggered bond configuration 
should be excluded from this list for \h\ > J 2 + J3 + ^ , since it cannot be the 
ground state due to the variational principle. 
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